AD-A062  727 


UNCLASSIFIED 


ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION  VICKSBURG  MISS 
FEASIBILITY  STUDY  OF  A NUMERICAL  TOW  MODEL. (U) 

SEP  78  T D ANKENY.  C J HUVAL » L L DAGGETT 
WES-MP-H-78-11  * 


I * I 


MISCELLANEOUS  PAPER  H-78-II 


FEASIBILITY  STUDY  OF  A 
NUMERICAL  TOW  MODEL 


Thomas  D.  Ankeny,  Carl  J.  Huval,  Larry  L.  Daggett 
Hydraulics  Laboratory 

U.  S.  Army  Engineer  Waterways  Experiment  Station 
P.  O.  Box  631,  Vicksburg,  Miss.  39ld0 


September  1978 
Final  Report 


Approved  For  Public  Rvlease,  Distribution  Unlimited 


■T  QUAXITT  PRAOTICAEUS, 
!l>  KDC  COflteAJMBD  A . 


Assistant  Secretary  of  the  Army  (R6>D) 
Department  of  the  Army 
Washington,  D.  C 20310 


Prepared  for 


Under  Project  No.  4AI6II0IA9ID 


Destroy  this  report  when  no  longer  needed.  Do  not  return 
it  to  the  originator. 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DDC  CONTAINED  A SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (W'.en  Data  Entered) 


READ  INSTRUCTIONS 


REPORT  DOCUMENTATION  PAGE 


BEFORE  COMPLETING  FORM 

3 RECIPIENT'S  CATALOG  number 


1.  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO. 


Miscellaneous  Paper 


4.  TITLE  (and  Subtitle) 


Final  repart. 


FEASIBILITY  STUDY  OF  A NUMERICAL  TOW  MODEL 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERfsJ 


7.  AUTHORf»J 


Thomas  D. ^Ankeny,  Carl  J.  Huval  Larry  L 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


AREA  4 WORK  UNIT  NUMBERS 

Program  Element  61101A 
Project  No.  i*Al6llfll_A2j 
TaskJ 1 Work  Unit  Wo. 

12.  REPORT  DATE 


U.  S.  Army  Engineer  Waterways  Experiment  Station 
Hydraulics  Laboratory 

Miss.  39180  ' 1/ 


P.  0.  Box  631,  Vicksburg 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Assistant  Secretary  of  the  Army  (R8tD)  •'"l  / \ 

Department  of  the  Army  l / ' 

Washington,  D.  C.  20310 

14.  MONITORING  AGENCY  NAME  & ADDRESS^//  different  from  Controlling  OKI ce) 


13.  NUMBER  OF  PAGES 


15.  SECURITY  CLASS,  (of  thla  report) 


Unclassified 

15a.  DECLASSI  F|  CATION/ DOWN  GRADING 
SCHEDULE 


16-  DISTRIBUTION  STATEMENT  (of  thla  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  In  Block  20,  If  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  (Continue  on  reverae  aide  If  neceaaary  and  Identify  by  block  number) 


Barges 

Feasibility  studies 
Mathematical  models 


Model  basins 
Towboat  s 
Tows  and  towing 


, TRACT  fConitoue  on  rarer  mo  aitba  If  neceeaary  mod  Identify  by  block  number) 


The  objective  of  the  study  described  in  this  report  was  to  explore  the 
feasibility  of  developing  a numerical  hydrodynamic  model  of  a typical  push 
towboat-barge  combination  for  use  in  engineering  planning  and  design  studies 
Such  a model  might  be  used  to  simulate  tow  movements  in  restricted  waterways 
in  critical  river  reaches  such  as  bends,  bridges,  and  near  navigation  locks 
and  dams  to  determine  the  adequacy  and/or  economic  efficiency  of  channel 
designs 


( Continued) 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


inclassified 

SECURITY  CLASSIFICATION  OF  THU  PAOEflWnn  DM*  Bnlmnd) 

20,  ABSTRACT  (Continued). 

A literature  search  was  conducted  and  it  was  found  that  no  information 
was  available  at  the  beginning  of  this  study  on  numerical  models  of  tow 
hydrodynamics  or  on  towing  tank  or  prototype  tests  of  the  maneuverability 
of  tows.  During  the  course  of  the  study  several  numerical  models  did  become 
available  and  further  research  in  this  area  was  initiated  by  the  U.  S.  Coast 
Guard  Research  and  Development  office.  This  work  is  described  herein. 

Since  no  information  was  available  for  the  development  or  testing  of  a 
numerical  tow  maneuvering  model,  a series  of  measurements  were  made  of  radio- 
controlled  scale  model  tows  used  in  physical  model  studies  at  the  Waterways 
Experiment  Station  as  they  executed  standard  ship  hydrodynamic  maneuvers. 

Data  from  these  measurements  and  computations  were  used  to  determine  estimates 
of  hydrodynamic  coefficients  of  a linear  model  of  ship  maneuverability. 

^Tt  was  found  that  it  is  feasible  to  predict  tow  maneuvers  using  a, 
numerical  model  of  tow  maneuverability;  however,  additional  measurements  of 
tow  response  characteristics  are  required  before  a model  that  would  be 
useful  in  engineering  studies  can  be  developed.  Comparisons  of  such  a model 
with  recent  measurements  of  full-scale  tow  maneuvers  could  demonstrate  the 
validity  of  this  model. 

\ 


Unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfUTiFn  Dmlm  E nfr»d) 


PREFACE 


The  investigation  reported  herein  was  conducted  under  Department 
of  the  Army  Project  No.  4ai61101A91D , In-House  Laboratory  Independent 
Research  (ILIR)  Program.  The  program  is  sponsored  by  the  Assistant 
Secretary  of  the  Army  (R8tD).  The  research  was  performed  at  the  U.  S, 
Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi. 

The  study  was  conducted  by  Messrs.  Carl  J.  Huval  and  Thomas  D. 
Ankeny  and  Dr.  Larry  L.  Daggett  of  the  Math  Modeling  Group  under  the 
general  supervision  of  Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics 
Laboratory  and  Mr.  M.  B.  Boyd,  Chief  of  the  Hydraulic  Analysis  Division. 
Mr.  M.  B.  Savage,  Instrumentation  Services  Division,  designed  the 
instrumentation  for  recording  the  tow  test  data  and  Mr.  L.  L.  Friar, 
Instrumentation  Services  Division,  assisted  in  recording  the  data  during 
the  tests. 

Directors  of  WES  during  the  conduct  of  this  investigation  and  the 
preparation  of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  John  L. 
Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 


CONTENTS 


Page 

PREFACE  1 

CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 

UNITS  OF  MEASUREMENT 3 

PART  I:  INTRODUCTION 1+ 

Background 1* 

Purpose  of  Study k 

PART  II:  REVIEW  OF  LITERATURE 7 

Research  Results  7 

Mathematical  Models  '7 

Other  Studies  in  Progress 10 

PART  III:  MATHEMATICAL  MODELS 13 

Linear  Model  13 

Nonlinear  Model  (Eda)  17 

PART  IV:  TOW  MODEL  TESTS 21 

Test  Facility  and  Conditions 21 

Model  Tow 22 

Test  Variables 22 

Data  Acquisition 27 

PART  V:  RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS  29 

Computation  and  Scale  Model  Results  29 

Study  Conclusions 32 

Recommendations  32 

REFERENCES 3** 

TABLES  1-11 

APPENDIX  A:  CALCULATION  OF  HYDRODYNAMIC  COEFFICIENTS  A1 

APPENDIX  B:  DEFINITION  OF  SYMBOLS  B1 


2 


t 


CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 


verted to  metric  (SI)  units  as  follows: 


Multiply 


degrees  (angle) 
feet 

feet  per  second 
inches 

knots  (international) 
pounds  (mass) 
square  feet 


-§5L 


0.0171*5329 

0.3048 

0.3048 

25.4 

0.5144444 

0.4535924 

0.09290304 


To  Obtain 


radians 

metres 

metres  per  second 
millimetres 
metres  per  second 
kilograms 
square  metres 
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FEASIBILITY  STUDY  OF  A NUMERICAL  TOW  MODEL 

PART  I;  INTRODUCTION 

Background 

1.  The  inland  waterway  system  is  composed  of  ports,  locks,  and 
the  connecting  channels  between  them.  As  the  system  has  grown  both  in 
traffic  density  and  miles  of  navigable  waterway,  the  need  to  have  a 
viable  means  of  designing  navigation  improvements  to  the  channels  has 
become  more  acute.  Proposed  channel  improvements  have  been  tested  by 
physical  hydraulic  models  and  model  towboats  with  attached  barge 
flotillas.  Figure  1 shows  a photograph  of  one  of  the  model  tows  being 
tested.  This  type  of  testing  and  analysis  has  been  useful  for  compara 
tive  channel  improvements  giving  a good  indication  of  tow  response, 
especially  in  high  current  regions.  The  reproducibility  of  tow 
navigability  results  is  difficult  due  to  the  model  piloting  variables 
that  affect  the  tests.  Physical  scale  modeling  is  also  very  expensive 
due  to  the  high  model  construction  cost. 

Purpose  of  Study 

2.  The  numerical  tow  model  project  was  initiated,  in  part,  to 
assess  the  benefits  that  would  result  from  channel  improvements 
(straightening,  widening,  or  lengthening  radii  of  bends)  over  a major 
segment  of  the  inland  waterway  system  instead  of  or  in  addition  to 
replacement  or  upgrading  of  lock  facilities.  It  is  presently  not 
feasible  to  build  a scale  model  to  study  the  impact  of  channel  improve 
ments  on  navigation  benefits.  Economic  analysis  of  channel  improve- 
ments requires  the  computation  of  corresponding  decrease  in  transit 
time.  The  relative  use  of  steering  and  flanking  rudders  and  towboat 
power  setting  also  influences  the  economy  of  tow  operation.  The 
purpose  of  this  study  was  to  explore  the  feasibility  of  developing  a 
mathematical  model  to  simulate  tow  behavior. 
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3.  This  report  describes  the  work  that  has  been  completed  on  the 
project  to  study  the  feasibility  of  developing  a numerical  hydrodynamic 
model  of  the  towboat-barge  combination.  The  report  covers  a brief 
assessment  of  relevant  literature  01  mathematical  ship  models  and 
the  efforts  to  develop  numerical  models  for  inland  tows.  A series  of 
tests  using  a scale  tow  model  are  described  and  the  resulting  data 
analysis  is  also  presented.  The  numerical  model  used  is  described  and 
the  results  of  typical  ship  and  tow  simulations  are  presented.  Methods 
of  calculating  the  required  tow  hydrodynamic  coefficients  are  discussed 
in  Appendix  A. 


PART  II:  REVIEW  OF  LITERATURE 


Research  Results 


h.  Some  of  the  results  of  research  on  slip  hydrodynamics  has  been 
presented  in  the  book  edited  by  Comstock.’*'  Ship  design  has  been 
improved  by  using  test  results  from  large  ship  model  towing  tanks  where 
data  on  hull  forms,  rudders,  propellors,  ship  stability,  and  ship 
maneuverability  are  systematically  developed.  The  mathematical 
modeling  of  ships  has  been  fairly  well  developed  and  reports  on  the 
results  of  such  research  are  numerous.  Most  ship  modeling  has  been 
concerned  with  ship  design  problems  in  open  water  navigation;  however, 
shallow  water  and  bank  effects  have  been  included  in  some  of  the  more 
recent  mathematical  models. 


Mathematical  Models 


5.  The  mathematical  models  developed  for  ships  are  based  on 
general  hydrodynamic  principles  and  are  equally  applicable  for  use  in 
a numerical  tow  model.  The  hydrodynamic  coefficients  that  are  re- 
quired in  the  numerical  model  are  expected  to  be  substantially  dif- 
ferent than  those  found  for  ships.  Some  examples  of  differences 
which  may  influence  the  numerical  model  are  given  below.  Towboat 
steering  systems  usually  require  a higher  level  of  maneuverability 
than  ships.  Combination  steering  and  flanker  rudders  with  Kort  noz- 
zles have  been  developed  for  towboats  involved  in  river  service. 
Typical  river  barge  flotilla  bows  are  much  blunter  than  ship  type 
forms  and  operate  at  much  shallower  drafts  than  ships.  The  large 
variability  of  barge  flotilla  configuration  is  another  factor  to 
consider  in  a numerical  model.  Towboat  propulsion  systems  are  also 
designed  considerably  different  from  ships  with  tunnel  hulls  and 
higher  power /load  ratios. 

6.  A review  of  the  literature  on  mathematical  models  of  towboat 
maneuverability  has  produced  very  few  studies  in  this  area;  at  the 
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beginning  of  the  study,  no  work  towards  developing  a mathematical  model 
of  tows  had  been  reported.  This  review  has  shown  that  towboat  design 
has  evolved  without  any  detailed  research  on  maneuverability,  such  as 
that,  with  ship  design.  This  may  be  due  to  the  fact  that  river  push 
towing  has  developed  relatively  recently  compared  to  ocean  shipping. 

There  appears  to  have  been  a lack  of  incentive  in  the  industry  to 
investigate  the  maneuvering  characteristics  of  barge  flotillas.  It 
was  learned  by  correspondence  and  personal  contacts  with  some  of  the 
representatives*  in  the  towboat  construction  industry  that  some  research 
work  has  been  accomplished  at  the  Netherlands  Ship  Model  Basin  to 
improve  towboat  propulsion  and  steering  efficiency.  Some  design  improve- 
ments resulting  from  these  studies  include  the  development  of  tunnel 
hulis,  Kort  nozzles,  large  engines  (up  to  10,000  hp),  large  towboats, 
and  integrated  tows.  Some  of  these  improvements  have  enhanced  the 


[.  • . . - maneuvering  characteristics  of  the  tows;  but  there  apparently  has  been 

little  effort  by  the  towing  industry  towards  systematic  study  of  tow 
maneuver ab i 1 i ty . 

2 

7.  Petrie's  towboat  and  barge  flotilla  model  became  available 
during  the  course  of  this  study.  This  model  appears  to  be  the  first 
to  address  tow  maneuverability  and  is  reviewed  in  some  detail  below. 

The  model  is  configured  to  predict  transit  times  through  a series  of 
maneuvers  for  various  barge,  towboat  and  bow  thruster  combinations. 

It  is  suited  for  assessing  the  influence  of  several  parameters  describ- 
ing the  physical  characteristics  of  the  tow,  such  as  its  length,  width, 
towboat  horsepower,  and  the  presence  or  absence  of  a bow  thruster.  The 
channel  through  which  the  tow  is  to  maneuver  can  be  described  by  a 
combination  of  straight  and  circular  segments,  and  a uniform  current 
velocity  can  be  specified.  The  model  does  not  simulate  all  of  the 
hydrodynamic  tow  characteristics  and  channel  depth  variations;  shallow 
water  and  bank  suction  effects  are  not  included.  The  mathematical 
model  integrates  the  three  force  and  moment  equations  in  the  horizontal 


* Dravo  Corporation,  Nashville  Bridge  and  Iron  Co.,  Netherlands  Ship 
Model  Basin,  personal  communication,  1977. 
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plane  to  determine  motion  and  position  of  the  flotilla.  The  hydrodynamic 
coefficients  used  in  the  model  can  be  either  based  on  towing  tank  test 
results  or  calculated  by  first  order  approbations  which  ignore  non- 
linear effects.  The  model  has  not  been  verified  against  laboratory 
or  field  data. 

8.  The  Petrie  model  was  procured  and  studied  for  possible 
adaptation  to  the  numerical  tow  model  development.  It  was  concluded 
that  this  model  could  be  of  some  use  in  river  planning  studies.  However, 
the  determination  of  hydrodynamic  force  and  moment  coefficients  would 
require  refinements.  Model  comparisons  and  validation  against  fielcl  or 
model  results  would  be  required.  A much  simpler  model  was  desired  for 
this  feasibility  study  of  tow  maneuvering  simulation.  However,  the 
model  represents  a useful  step  towards  the  eventual  goal  of  studying 

the  economic  impact  of  waterway  channel  improvement.  As  indicated  in 
Appendix  A,  the  hydrodynamic  coefficients  were  estimated  for  this  study 
using  Petrie's  calculation  method. 

9.  Recent  information  from  the  Coast  Guard  indicates  that  the 
Petrie  model  is  being  modified  in  the  following  areas: 

a.  Develop  more  generalized  channel  description  with  variable 
cross  currents, 

b.  Modify  propellor  and  rudder  model  to  allow  astern 
operation, 

£.  Improve  hydrodynamic  coefficients  estimations, 
d.  Incorporate  variable  wind  field  and  force  on  tow,  and 
e_.  Develop  a manual  pilot  control  capability. 

The  modified  model  software  will  be  available  to  the  Coast  Guard  at  the 
end  of  the  development. 

3 

10.  Eda's  ship  motion  simulation  model  also  became  available 

and  was  procured  during  the  course  of  this  study  for  possible  adaptation. 
The  model  formulation  is  presented  along  with  typical  tanker  maneuver- 
ing simulation  results.  Methods  to  generalize  the  model  to  include 
shallow  water,  bank  suction,  current  and  wind  effects  are  outlined  in 
the  report.  It  is  understood  that  these  efforts  are  presently  under 
way.  A suggested  set  of  hydrodynamic  coefficients  for  a typical 
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80,000  DWT  tanker  are  presented.  The  details  of  the  ship  numerical  model 

4 

are  described  by  Eda  , including  computer  code  listings.  Two  programs 
are  available;  one  program  simulates  standard  maneuvers  while  the  other 
simulates  typical  harbor  entering  maneuvers.  New  York  Harbor  is  used 
for  input  as  an  example  of  harbor  entrance  computations. 

11.  The  details  of  the  mathematical  ship  motion  model  developed 
by  Eda  were  studied.  It  was  concluded  that  the  model  for  the  standard 
maneuvers  (TURNCG)  could  be  adapted  to  the  towboat  numerical  modeling 
development.  Accordingly,  the  listed  computer  code  was  implemented  on 
the  Waterways  Experiment  Station  (WES)  timesharing  system.  The  model 
can  be  used  for  a towboat  as  well  as  ship  motion  simulation  of.  standard 
maneuvers  given  the  appropriate  hydrodynamic  coefficients. 

12.  The  program  NYHARB  for  harbor  entrance  modeling  is  longer 
and  somewhat  more  complex,  but  could  also  be  implemented  for  specific 
channel  waterway  studies.  The  model'’,  as  modified  by  the  Coast  Guard, 
can  simulate  the  behavior  of  a vessel  negotiating  a user  defined  harbor 
or  channel  system.  The  program  produces  a report  identifying  the 
periodic  location  and  status  of  the  simulated  vessel  over  time.  The 
model  incorporates  shallow  water  and  bank  suction  effects.  Because 
the  ship  must  negotiate  a harbor  system,  a steering  or  pilot  module 

is  included  to  provide  rudder  commands.  Elements  within  the  steering 
module  permit  the  inclusion  of  probabilistic  navigational  errors  if 
desired.  The  effects  of  wind  or  current  are  presently  not  incorporated 
in  the  program;  nor  is  squat  simulated.  Presently  the  hydrodynamic 
coefficients  are  available  only  for  an  80,000  DWT  tanker,  but  work  is 
under  way  to  provide  the  capability  to  simulate  a variety  of  types  and 
sizes  of  ships. 


Other  Studies  in  Progress 

13.  Recent  contact  with  the  Coast  Guard  also  indicates  a river 
tow  numerical  model  development  effort  under  way  at  Hydronautics , Inc. 
The  research  includes  the  following  tasks : 
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a.*  Modify  existing  ship  maneuvering  mathematical  models  to 
simulate  tow  motions, 

t>.  Planer  motion  mechanism  (PMM)  tests  to  determine  hydro- 
dynamic  coefficients  for  an  available  tow  model,  and 

£.  Simulation  study  of  river  tow  maneuvering  at  certain 
specific  bridge  sites. 

14.  The  model  as  presently  developed  is  for  deep  water  and  does 
not  simulate  shallow  water  or  bank  suction  conditions,  but  does  allow 
cross  velocity  shear  on  the  tow.  The  model  will  allow  both  ahead  and 
astern  thrusts  and  therefore  includes  flanking  maneuvers.  The  basic 
model  software  is  proprietary  and  will  not  be  available  under  the 
present  research  contract  with  the  Coast  Guard. 

15.  Due  to  the  lack  of  full  scale  tow  or  model  maneuverability 

data,  a cooperative  effort  was  initiated  for  WES  to  participate  through 

contract  in  a series  of  tow  performance  tests.  The  tow  performance  test 

program  was  expanded  to  include  tow  maneuvers  and  the  measurement  of 

parameters  that  would  provide  data  on  tow  maneuverability  characteristics. 

While  not  a specific  activity  of  the  tow  modeling  project,  it  is 

mentioned  briefly  here  to  show  how  the  WES  modeling  effort  led  to 

significant  advances  in  basic  tow  performance  information.  These  tow 

performance  data  can  be  used  to  validate  both  numerical  or  scale  tow 

models  for  use  by  WES  to  improve  navigation  simulation  validity.  These 

trials  were  conducted  on  the  Mississippi  River  at  Baton  Rouge, 

Louisiana,  in  November  1976.  A report  on  the  field  data  by  the 
7 

contractor  has  been  furnished  WES. 

16.  The  3360  hp  towboat  MV  Exxon  Memphis,  which  was  used  in 
these  tests,  is  powered  through  twin  screws  and  fitted  with  Kort  noz- 
zles. The  tow  flotilla  consisted  of  four  loaded  integrated  barges 
with  a total  length  of  ll60  ft  and  beam  of  5^  ft.*  Tests  were  con- 
ducted for  the  following  conditions: 

ji.  Steady  ahead,  straight  course  at  various  power  settings, 
up  river  and  down  river. 


* A table  for  converting  U.  S.  customary  units  of  measurement  to  metric 
(Si)  units  is  presented  on  page  3. 
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b.  Full  power  zig-zag,  up  river  and  down  river,  and 
c_.  Half  power  turns  up  river  and  down  river. 

Measurements  were  made  of  engine  speed  and  power,  shaft  speed,  and 
rudder  angle.  The  tow  position  and  attitude  were  obtained  from 
electronic  distance  measuring  equipment. 

17 . Additional  field  tests  were  conducted  in  November  1977  with 
two  additional  tows,  using  an  improved  data  collection  system  and 
involving  an  expanded  test  program.  It  is  expected  that  all  data  will 
be  published  and  distributed  for  use  by  the  industry  and  research 
laboratories . 
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derived  by  Abkowitz  and  the  details  are  presented  by  Mandel.  A 

method  of  solving  these  equations  for  ship  motion  with  a computer 

9 

program  has  been  presented  in  a David  Taylor  Model  Basin  report  for 
the  more  complete  nonlinear  model.  For  the  purpose  of  this  feasibility 
study,  the  linear  model  can  be  used  to  determine  whether  the  concepts 
used  for  ship  motion  simulations  can  be  used  for  towboat-barge  flotillas. 
The  equations  of  motion  are  formulated  for  a tow  with  rudders  for 
control  and  propellers  to  provide  thrust.  Tow  motion  is  translational 
in  the  x-y  plane  and  rotational  about  the  z-axis,  without  any  vertical 
motion.  It  has  been  found  that  these  three  equations  are  adequate  for 
maneuverability  simulation  of  large  surface  ships.  Figure  2 shows  a 
definition  sketch  and  the  coordinate  axis  about  the  tow  center  of 
gravity. 

19.  The  forces  and  moments  exerted  on  the  tow  are  composed  of: 

a..  Hydrodynamic  forces  and  moments,  due  to  the  motion  of 
the  tow, 

_b.  Rudder  and  propeller  forces,  and 

c_.  Disturbances  caused  by  wind,  waves,  current,  channel 
banks,  passing  tows,  etc. 

20.  For  this  study,  the  model  will  be  simplified  to  consider  the 
first  two  forces  only.  Disturbances  can  be  incorporated  after  further 
model  development.  The  model  under  consideration  will  be  formulated 
for  a tow  traveling  at  a constant  speed  in  deep,  slack  water  with  no 
cross  current  or  wind.  Hydrodynamic  and  rudder  forces  may  be  expressed 
in  a linear  form  by  expanding  the  equations  of  motion  in  a Taylor  series 
about  the  tow  center  of  gravity  and  considering  only  the  first  order 
terms.  Mandel^  gives  details  of  the  linear  model.  If  a tow  running 
with  only  small  deviations  from  a straight  course  at  constant  speed 
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is  considered,  two  linear  dimensional  equations  describing  the  sway  and 
yaw  motions  result:* 


(m  - Y.)  v - Y.r  - Y v - (Y  - inU)  r - Y.6  = EY.  (l) 

v r v r 6 1 


(I  - N. ) r - N.v  - N v - N r - N.6  = EN. 
z v v v r o 1 


where 


EY^,  EN^  are  external  forces  and  moments,  respectively, 
U is  the  speed  of  the  tow. 


Y.,  Y.  are  added  masses, 
v r 

N.,  N.  are  added  moments  of  inertia, 
r v 

Y , Y , N , N are  damping  coefficients,  and 
v r v r i-o 

Y.,  N_,  are  rudder  coefficients, 
o o 


Definitions  of  the  terms  are  given  in  Figure  2 and  Appendix  B.  In  a 
linear  model  the  axial  equation  (x  direction)  is  not  necessary  since  the 
tow  is  assumed  to  proceed  at  constant  speed  U so  that  U = . 

21.  The  dimensional  model  is  usually  made  dimensionless  and 
presented  as  follows  if  external  forces  and  moments  are  ignored: 

?t  iiit  i i i i 

(m  - Y.)  v'  - Y.r  - Y v - (Y  - m ) r = Yx6  (3) 

v r v r 6 

t If  If  If  II  I 

(I  - N . ) r - N.v  - H v -Nr  = N.6  (U) 

z r v v r o 

Definitions  of  the  dimensionless  hydrodynamic  coefficient  terms  are 
given  in  Table  1 and  Appendix  B. 

22.  The  linear  model  outlined  above  requires  the  specification  of 
at  least  nine  hydrodynamic  coefficients.  Theoretical  methods  can  be 


* For  convenience,  symbols  and  unusual  abbreviations  are  defined  in  the 
Notation  (Appendix  B). 
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used  to  obtain  some  of  the  coefficients,  but  towing  tank  tests  are 
required  for  reliable  tow  simulation.  A great  deal  of  testing  would  be 
necessary  to  measure  the  coefficients  to  produce  a fairly  complete 
linear  model  using  these  equations. 

23.  It  is  possible  to  simplify  the  mathematical  model  further  by 
using  the  Nomoto  equations.1^  These  equations  consider  the  hydrodynamic 
characteristics  of  the  tow  and  rudder (s)  to  be  characterized  by  two 
coefficients,  K and  T.  The  Nomoto  equations  are  much  simpler  than  the 
complete  linear  model  and  closed  form  solutions  have  been  obtained  for 
some  of  the  definitive  maneuvers.11  The  two  equations  for  yaw  and  sway 
are 

T»(t)  + «(t)  = K6(t)  (5) 

Tg  3 ( t ) + 3 ( t ) =•  K 6(t)  (6) 

where 

ip  is  the  heading  angle, 

$ is  the  drift  angle,  and 

T,  T.,  K,  K are  coefficients  based  on  the  hydrodynamics  of 

p p 

the  tow  and  rudder(s). 

The  drift  angle  b = - (in  radians)  for  a linear  system.  Dimensionless 

forms  of  these  equations  can  also  be  obtained.  Nomoto11  showed  how  the 

four  hydrodynamics  coefficients  K,  K„,  T,  T.  in  this  model  could  be  re- 
ts p 

lated  to  complex  relations  involving  the  linear  hydrodynamic  coefficients. 

2k.  The  attractiveness  of  the  simple  K-T  Nomoto  model  is  its 
simplicity  and  the  ability  to  obtain  the  four  coefficients  from  tow  field 
trials  or  free-running  scale  model  data.  Data  concerning  rudder  angle 
6,  drift  angle  $,  and  heading  angle  ip,  could  be  recorded  for  a particular 
maneuver  or  series  of  maneuvers.  By  reducing  the  data,  iji,  i|>,  and  6 can 
be  found  for  particular  values  of  6.  Using  data  from  several  tests  and 
the  Nomoto  equations,  a system  of  simultaneous  linear  equations  can  be 
solved  to  produce  the  coefficients  (K,  K,  T,  and  T0).  Once  the 

p p 
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coefficients  have  been  calculated,  the  solution  of  the  equations  for 
surge  and  sway  as  a function  of  rudder  settings  can  be  obtained  by  well 
known  finite  difference  time  stepping  procedures  if  initial  conditions 
are  known.  However,  to  produce  reliable  K-T  coefficients  would  require 
a rather  large  volume  of  model  data. 

Nonlinear  Model  (Eda) 

25.  The  complete  linear  model  and  the  simpler  Nomoto  model 
version  of  the  linear  model  are  applicable  only  at  small  rudder  and  yaw 
angles.  More  complete  nonlinear  models  have  been  formulated  and  used 
to  simulate  ship  motions  under  a variety  of  maneuvering  conditions. 

The  derivation  and  method  of  model  formulation  is  presented  in  detail 
by  Mandel\  Several  fairly  complete  nonlinear  models  are  known  to  have 
been  developed  using  square  or  cubic  polynomials  to  approximate  the 

3 4 

hydrodynamic  force  and  moment  functional  relations . Eda  ’ has 
presented  details  of  a model  for  ship  maneuvering  simulation  that  was 
recently  developed.  As  noted  above,  the  TURNCG  program  was  implemented 
on  WES  timesharing  system.  The  model  will  execute  either  a turning 
circle  or  z-maneuver  with  a pre-specified  sequence  of  rudder  angles. 
Figure  3 shows  numerical  ship  model  computational  results  for  an 
80,000  DWT  tanker. 

26.  The  mathematical  model  TURNCG  requires  several  kinds  of 
input  data  to  simulate  turning  circle  or  z-maneuvers  of  tows.  The 
program  allows  a choice  of  the  kinds  of  maneuvers  and  rudder  angle 
settings.  Table  2 gives  the  details  of  the  input  data  requirements. 

The  ship  hydrodynamic  coefficients  and  principal  dimensions  are  required 
as  input  data.  Initial  conditions  and  certain  control  parameters  to 
stop  the  computations  are  also  part  of  the  input  data. 

27.  The  computational  cycle  of  the  model  calculates  the 
hydrodynamic  moment  and  forces  at  each  time  step.  From  these  dimensional 
forces  and  moments,  the  linear  and  angular  accelerations  are  given  as 
follows : 


NOTE: 


LATERAL  DISTANCE.  yQ  = yQ.  L 

COMPUTATIONS  FOR  00,000  DWT  TANKER,  L = 763  FT. 
APPROACH  SPEED  = 16  KNOTS 


Figure  3.  Results  of  turning  circle  computations 


UDOT 

XT 

MX 

(7) 

TOOT 

YT 

MY 

(8) 

ROOT 

NT 

IZ 

(9) 

where  XT,  YT,  and  NT  are  the  sum  of  all  the  hydrodynamic  forces  and 
moments  acting  on  the  tow  and  MX,  MY,  and  IZ  are  the  tow  mass  plus  added 
mass  due  to  acceleration  or  tow  inertia  in  yaw  plus  added  moment  of 
inertia.  These  can  he  integrated  by  an  Euler  forward  difference 
technique  to  calculate  the  linear  and  angular  velocities  as  follows: 


U = U.  + UDOT*DELT 
x-1 


V = V.  + TOOT*DELT 
l-l 


R = R.  + RDOT*DELT 
l-l 


(10) 

(11) 

(12) 


The  tow  attitude  is  given  by  the  following  since  R = PSIDOT, 


PSI  = PSI.  . + R*DELT 
l-l 


The  distance  traveled  in  relation  to  the  fixed  axis  x , y by  the 

o o 

tow  can  be  calculated  from 


(13) 


X = X 
o 


/ . . N + UF*DELT 

ou-1) 


Yo  = Yo(i-l)  + VF*DELT 


(1U) 

(15) 


where 


UF  = U cos  PSI  - V sin  PSI  (l6) 

VF  = U sin  PSI  + V cos  PSI  (17) 


28.  The  output  of  the  model  gives  the  following  quantities  (see 
Table  3 for  an  example  of  the  model  output) 

TIME  - elapsed  time  from  beginning  of  maneuver,  sec, 

DELDEG  - rudder  angle,  in  degrees, 

PSIDEG  - heading  angle,  in  degrees, 

YP  - dimensionless  lateral  distance,  y^/L, 

XP  - dimensionless  longitudinal  distance,  xq/L, 

RP  - dimensionless  rate  of  turn,  r'  = rL/U  (note  that 
r is  in  radians  per  second  and  not  degrees  per 
second) , 

BETA  - drift  angle  from  tow  speed  vector  to  tow  axis,  in 
degrees , 

U/UO  - dimensionless  speed  as  a function  of  initial  straight 
course  speed  (speed  loss),  and 
TP  - dimensionless  time,  t'  = tU/L. 


The  output  file  could  be  easily  modified  to  include  ship  path  plots 
and  time  histories  of  tow  variables. 


PART  IV:  TOW  MODEL  TESTS 

29.  The  literature  study  and  personal  contacts  with  various 
companies  and  institutions  showed  that  no  data  on  towboat  maneuverability 
was  available.  Thus,  it  was  decided  to  conduct  a series  of  tests  in 

the  laboratory  using  a model  towboat  barge  combination  in  a deep  water 
test  basin  for  slack  water  condition.  Slack  water  conditions  are 
different  from  river  operations  which  have  strong  to  moderate  current 
patterns;  however,  these  conditions  were  considered  adequate  since  the 
purpose  of  the  present  study  was  to  determine  the  applicability  of 
hydrodynamic  models  of  ships  to  shallow  draft  tows.  The  model  test 
results,  while  limited,  would  allow  the  development  of  a first  stage 
numerical  towboat  model  to  be  refined  as  data  become  available. 

Test  Facility  and  Conditions 

30.  This  section  describes  the  test  facility,  test  conditions, 
model  tow,  variables  instrumented,  and  data  acquisition.  Tow  tests 
were  conducted  with  two  of  the  ship  type  standard  maneuvers , the  circle 
test  being  the  primary  maneuver.  These  tests  can  be  used  to  determine 
the  performance  of  the  tow  at  various  rudder  settings.  Testing  was 
conducted  in  a large  tidal  inlet  physical  model  facility  at  WES  during 
February  and  March  1977-  The  overall  dimensions  of  the  testing  basin 
were  60  ft  x 128  ft.  The  actual  area  used  varied  from  test  to  test 
depending  on  the  particular  maneuver  being  simulated.  The  water  level 
was  maintained  so  that  the  depth  during  testing  was  7-1/2  to  8 in.  A 
reference  grid  system,  composed  of  numbers  (y-axis)  and  letters  (x-axis), 
was  painted  on  the  horizontal  concrete  bottom  with  an  interval  of  2 ft. 
The  basin  was  located  in  an  inclosed  building  where  no  wind  was  present. 
An  observation  area  for  the  video  tape  camera  operator  was  provided  by 
cat  walks  above  the  model  area. 
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Model  Tow 


31.  The  model  tow  consisted  of  a towboat  and  flotilla  of  barges 
as  used  in  physical  model  studies  and  is  shown  in  Figure  1.  Dimensions 
of  the  towboat  and  barge  flotilla  are  given  on  Table  4.  The  towboat 
was  constructed  of  plastic  having  two  propellers  fitted  with  Kort 
nozzles,  with  a steering  rudder  aft  and  a flanking  rudder  forward  of 
each  propeller.  Figure  4 presents  a drawing  of  the  model  towboat.  The 
propellers  were  driven  by  an  electric  motor  whose  shaft  was  geared 
directly  to  each  propeller  shaft.  The  motor  could  be  controlled  at  a 
constant  speed  either  in  the  forward  or  reverse  direction  by  a remote 
control  unit.  The  rudders  operated  as  a unit  (steering  or  flanking) 
and  were  also  actuated  by  the  remote  control  unit.  The  flotilla  was 
modeled  as  a single  unit  integrated  barge.  The  dimensions  are  given 
on  Table  4 and  a definition  sketch  presented  in  Figure  5.  The  barge 
was  raked  at  the  bow  with  a smooth,  flat  bottom  and  straight  sides  and 
stern,  typical  of  modern  integrated  barges.  The  barge  was  constructed 
of  sheet  metal  and  held  the  power  supplies  for  the  boat  and  the 
instrumentation  as  well  as  the  data  transmitter.  Additional  ballast  was 
provided  by  bricks  that  were  positioned  to  balance  the  barge  both 
longitudinally  and  transversly.  This  towboat-barge  combination  has  been 
used  by  the  Waterways  Division  of  the  Hydraulics  Laboratory  to  simulate 
a 600  ft  x 105  ft  tow  pushed  by  a 5000  hp  towboat  (approximate  scale 
1:70)  for  studies  at  lock  approaches  in  bendways,  and  other  navigation 
problems.  During  testing  only  the  forward  drive  and  the  steering 
rudders  were  used.  The  flanking  rudders  remained  in  a neutral  position. 

Test  Variables 

32.  The  data  items  recorded  both  on  magnetic  tape  and  strip 
charts  were:  angular  position  of  steering  and  flanking  rudder(s),  speed 
of  the  motor  in  rpm,  and  the  motor  current  and  voltage.  Rudder  angles 
were  obtained  by  recording  voltages  from  constant  rotation  potentiometers 
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towboat  characteristics 


attached  to  the  steering  mechanisms.  Speed  of  the  motor  was  from  the 
tachometer  built  into  the  boat  by  the  manufacturer. 

33.  The  testing  program  was  composed  of  a series  of  straight  line 
runs  to  determine  the  velocity  of  the  tow  for  various  motor  settings; 

a series  of  constant  rudder  angles  and  speed  settings  which  produced 
the  circle  test,  and  a zig-zag  run.  More  details  of  each  segment  of 
the  test  program  are  given  below  and  are  summarized  in  Table  5- 

34.  The  straight  line  runs  to  calibrate  the  speed  of  the  tow 
were  made  by  selecting  a throttle  setting,  lining  up  the  tow,  starting 
forward  motion  and  then  timing  the  transit  between  markers  after  the 
tow  had  reached  constant  speed.  During  the  test  only  the  steering 
rudders  were  used  to  maintain  the  straight  course.  Four  throttle 
settings  were  chosen  and  a complete  test  was  composed  of  both  a North- 
South  run  and  a South-North  run.  Data  from  these  tests  are  presented 
in  Table  6 and  plotted  on  Figure  6. 

35.  Combinations  of  three  throttle  settings,  six  steering  rudder 
settings,  (three  angle  settings  both  starboard  and  port)  and  three  trials 
per  set  of  conditions  were  run  to  produce  a large  number  of  circle  test 
data.  The  steering  rudder  was  set  using  a stop-plate  to  maintain  a 
constant  rudder  angle  during  the  entire  circle  test.  The  throttle 
setting  was  constant  during  the  circle  tests.  The  circle  tests  were 
conducted  by  propelling  the  tow  on  a straight  line  course,  steady  ahead 
with  a sudden  command  changing  the  rudder  angle  to  the  desired  setting. 
Large  rudder  settings  produced  small  diameter  circles  and  a complete 
circle  was  able  to  be  completed  within  the  test  basin.  For  some  of 

the  smaller  rudder  settings,  the  circle  was  so  large  that  only  a portion 
(approximately  1/2  to  3/4  of  a circle)  could  be  completed  due  to  the 
limited  physical  size  of  the  facility. 

36.  The  zig-zag  maneuver  consisted  of  initiating  tow  motion  along 
a straight  course,  changing  the  steering  rudder  setting  to  approximately 
10  deg  in  one  direction,  waiting  for  the  tow  to  attain  a 10-deg  attitude 
change  and  then  changing  the  rudder  angle  in  the  opposite  direction. 
Repetition  of  this  procedure  constituted  a zig-zag  maneuver.  The  change 
in  attitude  was  by  visual  observation  only. 
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Data  Acquisition 


37-  All  tests  were  recorded  on  video  tape  to  obtaii  position  and 
attitude  of  the  tow,  and  on  magnetic  tape  and  strip  chart  paper  to  record 
steering  rudder  angle,  flanking  rudder  propeller  angle,  aid  motor  voltage 
and  current.  Two  specific  tests  were  analyzed  in  detail  for  this  report. 
The  remaining  test  data  are  available  for  further  analysis  , if  desirable. 

38.  The  video  tape  was  taken  using  a hand-held  camera  approximately 
18  ft  above  the  water  surface.  A time  code  generation  de\ ice  was 
incorporated  with  the  taping  so  that  the  video  could  be  ccrrelated  with 
the  other  recorded  data.  During  the  testing,  it  was  fount  that  the 
video  time  code  generator  was  not  in  perfect  synchronization  with  the 
code  on  the  magnetic  and  paper  tapes.  The  video  indication  of  t im< 
was  about  lit  sec  slow,  i.e.,  an  event  recorded  in  slow  code  on  the 
strip  chart  recorder  at  1:00.00  would  be  recorded  on  video  with  a ’ im< 
indication  of  l:00.lU. 

39*  The  data  recorded  on  the  strip  chart  from  the  tow  was  a 
function  of  a voltage  measurement.  Zeroes  and  calibration  steps  pre  ■.  ded 
each  series  of  tests.  The  recorded  voltages  were  translated  t„  numerical 
values  of  the  pertinent  tow  variable  with  the  use  of  calibration  curves. 
The  chart  data  were  analyzed  from  the  strip  chart  record  with  samples 
every  10  sec.  Visual  inspection  of  the  chart  records  revealed  nearly 
constant  values  for  the  tow  variables  thus  more  frequent  measurements 
were  not  required.  The  reduced  data  for  the  two  tests  that  were  analysed 
are  presented  in  Tables  7 and  8. 

40.  Analysis  of  the  video  tape  to  obtain  the  pos:tion  and 
orientation  of  the  tow  during  the  test  maneuvers  was  rather  difficult 
and  time  consuming.  This  was  due  to  the  lack  of  adequate  contrast  on 
the  video  tape  and  difficulty  in  stopping  the  tape  at  a desired  time 
for  reading  the  tow  position.  Extraction  of  the  position  of  the  tow's 
bow,  stern  and  center  of  gravity  was  obtained  from  still  photographs 
taken  at  various  times  during  the  test  maneuvers  as  recorded  on  the 
video  tape.  The  positions  were  measured  with  respect  to  the  grid  that 
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had  been  created  on  the  test  basin  floor.  These  grid  locations  were 
translated  into  x-y  coordinates  and  calculations  of  , 4*,  U,  and  3 
were  based  on  these  coordinates  and  rates  of  change  in  the  coordinates 
using  simple  geometric  relationships.  Tables  9 &nd  10  present  the 
reduced  data  from  the  video  tape  recordings. 
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PART  V:  RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 


Computation  and  Scale  Model  Results 


Ul.  The  data  from  the  reduced  scale  model  test  No.  IT  are  shown 
on  Figure  7 compared  with  the  numerical  tow  model  computations.  As  ex- 
plained in  Appendix  A,  the  rudder  force  and  moments  were  corrected  on  the 
basis  of  the  observed  model  data.  In  addition,  as  shown  in  Appendix  A, 
only  linear  tow  coefficients  could  be  calculated;  thus,  the  computations 
result  from  a linear  model.  The  coefficients  used  to  obtain  the  results 
shown  in  Figure  7 are  those  resulting  after  all  adjustments  an;  are  the 

last  values  listed  for  each  coefficient  in  Table  1.  The  coefficients 

L 

were  used  in  Eda's  TURNCG  model  to  calculate  the  tow's  turning  circle. 

1+2.  The  linear  mathematical  model  can  be  used  to  find  several 

important  stability  and  steering  qualities  of  ships  and  tows.  in 

particular,  the  relation  between  a turning  circle  radius  and  the 

■] 

hydrodynamic  coefficients  have  been  derived'1'  and  are  given  below; 


R 


L 

6 


Y N.  - N I. 
v 6 vo 


(18) 


i i 

Usually,  N.  and  Y.  are  given  in  moment  and  force  per  radian  so  that 
o o 

6 has  to  be  given  also  in  radians.  The  dimensionless  rate  of  turn 


Ij  _ rL  _ jjdj 
R " U = U 


(19) 


The  values  of  r = ip  are  also  used  in  radians. 

1+3.  The  equation  for  the  radius  of  the  turning  circle  was  used 

» » 

to  obtain  improved  values  of  L and  Yt  to  fit  the  observed  tow  model 

o o 

turning  circles.  If  the  tabulated  values  of  the  tow  hydrodynamic 
coefficients  are  used  in  the  above  equation,  the  radius  and  other 
wide-turn  test  values  are  found  to  be 

R = 11+.39  ft,  D = 28.8  ft,  R/L  = 1.81,  r'  = 0.55f» 
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It  is  also  possible  to  calculate  the  rate  of  turn  using  the  yaw  angle 
as  a function  of  time  to  obtain  i|;.  The  tow  yaw  data  can  be  used  to 
obtain  the  following  angular  rate  of  turns: 


Test 

n ^ 

U 

1 

R 

No. 

( /sec) 

( ft/ sec ) 

r 

(ft) 

R/L 

ll* 

3.509 

1.18 

0.1*13 

19.2 

2.1*2 

IT 

1.835 

0.60 

0.1*25 

18.7 

2.35 

These  results  indicate  that  the  tow  speed  has  only  a minor  effect  on , 
rate  of  turn.  The  rate  of  turn  from  the  yaw  rate  data  indicates  a 
smaller  rate  of  turn  and  longer  turning  circles  by  about  30  percent 
compared  to  the  geometrically  determined  rate  of  turn.  This  is 
probably  due  to  speed  loss  effects  during  the  turning  circle  which  are 
not  properly  simulated  in  the  linear  model. 

2.2 

1*1*.  Hydrodynamic  coefficients  have  been  published  by  Eda  for  a 
barge  model  flotilla  with  a push  tow.  The  configuration  was  quite 
different  than  the  tow  scale  model  used  to  obtain  the  circle  test  data 
presented  in  this  report.  Nevertheless,  the  comparison  may  be  worth 
noting  and  is  presented  in  Table  11.  Also  shown  are  preliminary  towing 
tank  results  obtained  with  a typical  river  towboat  and  barge  flotilla. 1 
Also  shown  are  computations  of  turning  radius  and  rate  of  turn  for  a 
20°  rudder  angle.  Additional  research  effort  will  be  required  to 
properly  interpret  these  computations. 

( 

1*5.  One  of  the  Exxon  full  scale  tow  tests  previously  mentioned 
included  a steady  upriver  turn  with  the  following  characteristics: 

^ = -0.196°/sec, 

6 = +ll*.9°, 

L = ll60  ft,  and 

U I 13.1*  ft/sec  (speed  of  approach  through 
the  water ) . 

From  these  data  the  rate  of  turn  can  be  calcuated  as 


radians 


I 

i 

t 


^ = 0.296 


If  the  rate  of  turn  versus  rudder  angle  is  approximately  linear,  this 

* o 

would  give  r = 0.397  at  20  . 

h6.  The  comparative  tow  maneuverability  data  shown  above  are 
presented  to  show  the  preliminary  nature  of  the  numerical  model.  The 
need  for  more  definitive  information  on  tow  hydrodynamic  coefficients 
is  obvious.  However,  the  comparison  of  results  shown  on  Figure  7 
indicates  a reasonably  good  fit  of  the  scale  model  data  is  possible. 
Additional  computations  and  comparisons  to  other  model  tests  would  be 
necessary  to  more  completely  validate  the  numerical  model.  Nevertheless, 
the  results  demonstrate  the  feasibility  of  such  computations. 


Study  Conclusions 

1*7 . The  study  has  shown  that  there  were  no  known  reliable 
numerical  models  specifically  applicable  to  tow  motion  simulation  at 
the  time  of  this  study.  Mathematical  models  that  are  well  developed 
for  ship  maneuverability  can  be  adapted  to  tows , but  require  several 
tow  hydrodynamic  coefficients.  Towing  tank  tests  and  numerical  model 
development  are  being  conducted  under  Coast  Guard  sponsorship.  These 
measurements  and  the  hydrodynamic  coefficients  are  only  now  becoming 
available.  Standard  maneuver  tests  with  free  running  models  (such  as 
those  conducted  in  this  study)  can  be  valuable  in  validating  a devel- 
oped numerical  tow  model.  Based  on  the  results  of  this  exploratory 
model  development,  it  is  concluded  that  a useful  numerical  tow  model 
for  tow  motion  simulation  is  feasible. 

Recommendations 


U8.  The  ongoing  research  at  several  research  laboratories  should 
be  monitored  for  hydrodynamic  coefficient  test  results  and  improvements 
in  numerical  tow  modeling.  When  the  data  and  refined  models  become 
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available,  additional  development  work  should  be  undertaken  to  improve 
the  capability  for  tow  modeling.  It  is  recommended  that  numerical  tow 
models  be  validated  against  full-scale  trial  and  free-running  scale 
model  data. 


r 
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TABLE  1 


- 

p •' 


SUNMARY  OF  HYDRODYNAMIC  COEFFICIENTS 
Symbol Definition  Magnitude (xlO  3) Remarks 


Nv 


N' 

V 

0.0961 

Yaw  moment  due  to  sway 

V 

1/2  pL3U 

N’ 

r 

N 

*Nr 

A 

-0.5899 

Yaw  damping 

1/2  pL  U 

N I 

- I' 

r - z 

r 

-0.0487 

Added  yaw  moment  of  indrtia 

z 

1/2  pL3 

(-0. 3608) * 
(-0.2291)** 

N ' 

<5 

N6 

1/2  pL3U2 

-0.5695 

(-0.98)*** 

Rudder  moment 

Y 

Y 1 

V 

2 

-15.5325 

Sway  damping 

V 

1/2  pL  U 

Y 

Y ’ 
r 

r 

1/2  pL3U 

Y . - m 

V 

7 

0.0961 

Sway  force  due  to  yaw 

- m ' 

-0.8027 

Added  mass  in  sway 

1/2  pL0 

(-7.127)* 

(-4.3662)** 

Y. 

Y6 

5 

1/2  pL2U2 

0.9520 

(1.63)*** 

Rudder  force 

X' 

X6 

*>  2 

-C . 9520 

Rudder  drag 

6 

1/2  pL^TJ 

(-1.63)*** 

I1 

I 

z 

r 

0.1804 

Moment  of  inertia 

Z 

1/2  pL 

m ’ 

m 

7 

3.5635 

Mass  of  tow 

1/2  PL 

* Estimated  using  21'  and  2m'  , respectively.  I 

**  Calculated  assuming  Petrie's  formulas  calculated  Nj.  and  Y(,  , 
rather  than  N.'  - I'  , and  Y!  - m'  , (e.g.,  N!  = -0.0487  - 0.1804 

rz  V ° r 

= -0.2291). 

***  Corrected  on  the  basis  of  measured  model  turning  diameter 

(1.72  x N'  Y'  and  Xl). 

0 0 0 
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INPUT  DATA  INFORMATION 
Numerical  Tow  Model 


FORTRAN 


Symbol 

Def ini t i on 

SL 

L 

Total  Tow  length,  ft 

UI 

U. 

l 

Initial  speed,  ft/sec 

PRDEL 

At 

Time  increment  for  printing  trajectory  data,  sei 

FINTIM 

t 

Time  limit  to  finish  a computer  run,  sec 

FINPSI 

Heading  angle  limit  to  finish  a run,  degrees 

RI.B 

L/B 

Length-beam  ratio 

RBH 

B/H 

Beam-draft  ratio 

CB 

Cb 

Block  coefficient  = V/LBH  = A/yLBII 

CP 

C 

p 

Prismatic  coefficient  = V/A^L,  A^  = area  amidsinj  septic'. 

AR 

ar/lh 

Rudder  area  ratio 

RHDW 

H/Dw 

Draft/water-depth  (It  is  zero  in  deep  water. 1 

RBW 

B/W 

Beam/channel-width  (It  is  zero  in  open  water.) 

A (N) 

N=1 ,15 

Hydrodynamic  moment  coefficients  in  yaw  equation 

B(N) 

N-1,13 

Hydrodynamic  force  coefficients  in  sway  equation 

C(N) 

N=1 ,8 

Hydrodynamic  force  coefficients  in  axial  or  surge  equation 

NIDF.L 

-- 

Number  of  rudder  angles 

RUDI  (N) 

S 

Rudder  angle  in  degrees,  N=1 , NIDFI, 

EXTNO 

-- 

Execute  number  (1  and  5 for  circle  turn-  and  z-maneuver, 
respectively) 

RUDZ 

Heading  angle  in  degrees  which  initiates  rudder  execution 

to  the  opposite  direction  during  z-maneuver 


TABLE 
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TIME 

0. 

20.  00 

40.00 
60.  00 

80.00 

1 00. 00 

120. 00 

140. 00 

160. 00 

180.00 

200. 00 
220.:  00 

240. 00 

260. 00 

280. 00 

300. 00 

320. 00 

340. 00 

360. 00 

380. 00 

400. 00 

420.00  " 

440. 00 

460. 00 

480. 00 

500. 00 

520. 00 

540. 00 

560. 00 

580. 00 

600. 00 

620. 00 
^40. 00 

660. 00 

680. 00 

700. 00 

720. 00 

740. 00 

760. 00 


Example  of  Numerical  Model  Output  for  Turning  Circle 

80,000  DWT  Tanker,  L = 763  ft 
Approach  Speed  = 16  Knots 


DELDEG 

PS I DEG 

YP 

XP 

RP 

BETA 

LMJ0 

TP 

0. 

0. 

0. 

0. 

0. 

0. 

0.46 

0. 

-15.00 

T.33 

0.  00 

0.  66 

0.11 

1'.  49 

1.00 

0.66 

-15. 00 

7Y54 

0.  02 

1.32 

0.  19 

3.74 

0.99 

1.32 

-15. 00 

15.74 

0.  10 

1.97 

0.25 

5.91  • 

0.97 

1.97 

-15V00 

25.51 

0.25 

2.53 

0.29 

7.74 

0.94 

2.60 

—151 0 0 

36.04 

0.48 

3.  15 

0.31 

9.  15 

0.90 

3.21 

-t5i.'00 

46.85 

0.78 

3.64 

0.33 

0.21 

0.87 

3.30 

-15.00 

57..  68 

1.16 

4.  06 

0.34 

1. 03 

0.33 

4.36 

-15. 00 

63.  44 

1.58 

4.40 

0.35 

1.66 

0.79 

4.90 

-15. 00 

79.  07 

2.03 

4.64 

0. 36 

2.  17 

0.76 

5.42 

-15. 00 

89.55 

2.50 

4.80 

0.37 

2.59 

0.73 

5.91 

-15. 00 

99.88 

2.98 

4.86 

0.33 

2.92 

0.71 

6.39 

-15. 00 

110.07 

3.44 

4.85 

0.39 

3.  19 

0.69 

6.36 

-15. 00 

120. 11 

3.88 

4.76 

' 0.39 

3.41 

0.67 

7.31 

-15. 00 

130. 04 

4.29 

4.60 

0.4  0 

3.59 

0. 66 

7.75 

-15.00 

139.35 

4.66 

4.37 

0.40 

3.73 

0.64 

8.  18 

-15.00 

149.56 

4.98 

4.  10 

0.40 

3.84 

0.63 

8.60 

-15. 00 

159. 19 

5.24 

3.78 

0.40 

3.93 

0.62 

9.  02 

-15. 00 

168.75 

5.45 

3.42 

0.41 

3.99 

0.62 

9.43 

-15. 00 

173.25 

5.59 

3.  04 

0.41 

4.  05 

0.61 

9.34 

-15. 00 

137.69 

5.67 

2.65 

0.41 

4.  09 

0.60 

10.24 

-15. 00 

197. 09 

5.68 

2.25 

0.41 

4.  12 

0.60 

1 0.64 

-15. 00 

206.45 

5.63 

1.86 

0.41 

4.  14 

0.60 

11. 04 

-15. 00 

215.78 

5.52 

1.48 

0.41 

4.  16 

0 . 59 

11.43 

-15. 00 

225. 09 

5.35 

1.  13 

0.41 

4.  17 

0.59 

11.82 

-15. 00 

234.38 

5.  12 

0.81 

0.42 

4.  18 

0.59 

12.21 

-15. 00 

243.65 

4.85 

0.54 

0.42 

4.  19 

0.59 

12.6  0 

-15.00 

252.90 

4.53 

0.31 

0.42 

4.20 

0.53 

12.99 

-15. 00 

262. 15 

4.  19 

0.  13 

0.42 

4.20 

0.58 

13.33 

-15. 00 

271.38 

3.82 

0.  02 

0.42 

4.21 

0.58 

13.77 

-15. 00 

230.61 

3.44 

0.  04 

0.4? 

4.21 

0.58 

14.  15 

-15. 00 

289.83 

3.  05 

0.  08 

0.42 

4.21 

0.58 

14.54 

-15. 00 

299. 04 

2 . 67 

0.  03 

0.42 

4.22 

0.58 

14.92 

-15. 00 

303.25 

2.31 

0.  16 

0.42 

4.22 

0 . 53 

15.31 

-15. 00 

317.46 

1.97 

0 . 34 

0.42 

4.22 

0.53 

15.69 

-15. 00 

326. 66 

1.67 

0.  58 

0.42 

4.22 

0.58 

16.  03 

-15. 00 

335.87 

1.40 

0.86 

0.42 

4.22 

0.53 

16.  46 

-15. 00 

345. 07 

1.19 

1.18 

0.42 

4.22 

0.58 

16.85 

-15. 00 

354.26 

1 . 03 

1.53 

0.42 

4 . 22 

0.58 

17.23 

TABLE  4 

PRINCIPAL  PARTICULARS 
Model  Towboat  and  Flotilla 


(Approximate  Scale  1 

: 70) 

Dimension- 

less 

ft 

lbs 

sq  ft 

Tow  Length 

5.583 

Towbfcat  Langth 

2.384 

Total  Length 

7.967 

Tow  Beam 

1.500 

Towboat  Beam 

0.5208 

Tow  Draft 

0.1042 

Towboat  Draft 

0.09375 

Tow  Center  of  Gravity  (C.G.)  Aft 

2.813 

Towboat  C.G.  Aft 

1.271 

Flotilla  and  Towboat  C.G.  Aft 

Location  Tow  plus  boat  C.G. /Total  Length 

0.4011 

2.391 

Tow  Displacement 

50.3 

Towboat  Displacement 

5.8 

Total  Displacement 

56.1 

Tow  Block  Coefficient 

Towboat  Block  Coefficient 

Total  Area  Steering  Rudders  (Twin) 

Total  Area  Flanking  Rudders  (Twin) 

Steering  Rudder  Area  Coefficient  (^~) 

0.9240 

0.8016 

0.03847 

0.03201 

20.909 

Steering  Rudder  Aspect  Ratio  (Span/chord) 

0.S620 

Propeller  Diameter  (Twin  Screws, 

0.1083 

Three-bladed , tips  cut  off) 

Kort  Nozzle  Inside  Diameter  (Twin  Screws) 
Propeller  Pitch/Diameter 

1.814 

0.1160 
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TABLE  7 

Circle  Test  No.  14,  Speed  Settinp  = 14 
Reduced  Strip  Chart  Data 


Motor 

Motor 

Current 

Rudder  Angles 

Time 

(min: sec) 

Shaft 

RPM 

Volts 

HP 

xlO-3 

Steering 

n 

Flanker 

(°) 

3214 

3.92 

1.16 

6.098 

-19.0 

-4.3 

3148 

3.88 

1.12 

5.828 

-19.8 

-7.1 

00:40 

3247 

4.00 

1.18 

6.330 

-19.8 

-8.0 

00:50 

3247 

4.00 

1.16 

6.222 

-19.0 

-8.0 

01 :00 

3247 

3.96 

1.16 

6.160 

-19.8 

-7.1 

01:10 

3247 

3.96 

1.16 

6.160 

-19.8 

-7.1 

01:20 

3247 

3.96 

1.16 

6.160 

-19.8 

-8.0 

01:30 

3247 

3.96 

1.14 

6.054 

-20.7 

-8.0 

01 :40 

3247 

3.96 

1.14 

6.054 

-19.8 

-7.1 

01 :50 

3247 

4.00 

1.14 

6.115 

-19.8 

‘-5.3 

02:00 

3247 

4.00 

1.14 

6.115 

-20.7 

-7.1 

02:10 

3281 

4.00 

1.16 

6.222 

-20.7 

-7.1 

02:20 

3281 

4.00 

1.16 

6.222 

-20.7 

-3.3 

02:30 

3314 

4.04 

1.18 

6.393 

-19.8 

-5.3 

02:40 

3281 

4.00 

1.16 

6.222 

-19.8 

-6.2 

02:50 

3281 

4.04 

1.16 

6.285 

-20.7 

-7.1 

03:00 

3314 

4.04 

1.16 

6.285  . 

-19.8 

-7.1 

03:10 

3314 

4.04 

1.16 

6.285 

-20.7 

-7.1 

03:20 

3281 

4.04 

1.14 

6.176 

-19.8 

-8.0 

03:30 

3247 

4.00 

1.16 

6.222 

-19.8 

-7.1 

03:40 

3247 

3.92 

1:14 

5.993 

-19.8 

-8.0 

03:50 

3247 

3.96 

1.16 

6.160 

-20.7 

-7.1 

r~ 


1 


TABLE  8 

Circle  Test  No.  17,  Speed  Setting  = 8 
Reduced  Strip  Chart  Data 


Time  Shaft 

(min  :sec) RPM 


Motor  Motor 

Volts  Current 


Rudder  Angles 

Steering  Flanker 

(°) n 


39:30 

1624 

2.16 

0.58 

1.680 

3.5 

-5.3 

39:40 

1690 

2.12 

0.62 

1.763 

-20.7 

-7.2 

39:50 

1690 

2.16 

0.60 

1.738 

-20.7 

-7.2 

40:00 

1690 

2.16 

0.62 

1.796 

-20.7 

-7.2 

40:10 

1690 

2.16 

0.64 

1.854 

-20.7 

-7.2 

40:20 

1690 

2.20 

0.64 

1.888 

-19.8 

-7.2 

40:30 

1690 

2.16 

0.64 

1.854 

-19.8 

-5.3 

40:40 

1723 

2.20 

0.64 

1.888 

-19.8 

-6.3 

40:50 

1723 

2.20 

0.62 

1.829 

-19.8 

-7.2 

41:00 

1723 

2.16 

0.62 

1.796 

-19.8 

-7.2 

41:10 

1723 

2.16 

0.64 

1.854 

-19.8 

-5.3 

41:20 

1756 

2.24 

0.64 

1.922 

-19.8 

-7.2 

41:40 

1723 

2.16 

0.60 

1.738 

-20.7 

-6.3 

42:00 

1723 

2.16 

0.62 

1.796 

-19.8 

-2.6 

42:30 

1657 

2.16 

0.60 

1.738 

-20.7 

-6.3 

43:00 

1657 

2.16 

0.60 

1.738 

-20.7 

-6.3 

43:30 

1690 

2.16 

0.58 

1.680 

-20.7 

-6.3 

44:00 

1657 

2.08 

0.58 

1.618 

-19.8 

-5.3 

C4  tO 


TABLE  9 


il 


Circle  Test  No.  14,  Speed  Setting  = 14 
Tow  Position  from  Video  Tapes 

Position  of 

Center  of  Gravity  .. 

Time  At  ft  UCG 


(Min:  sec)  (Sec)  Y0  Xp  \|»*  ft/sec  B* 


00 

00 

28 

30 

2 

6.00 

4.00 

32.10 

0.0 

0.0. 

00 

34* 

4 

4.15 

37.17 

4.3 

1.27 

2.6 

00 

37 

3 

5.30 

40.84 

23.6 

1.28 

6.2 

00 

40 

3 

4.10 

44.36 

15.0 

1.24 

-3.9 

00 

43 

3 

4.40 

45.82 

27.2 

.50 

15.6 

00 

44 

1 

6.34 

50.40 

39.3 

4.97 

16.3 

00 

47 

3 

6.10 

50.88 

37.5 

.18 

11.0 

00 

58 

11 

14.65 

56.00 

90.0 

.91 

30.9 

01 

00 

2 

16.58 

58.21 

88.4 

1.47 

47.3 

01 

12 

12 

27.20 

56.29 

139.3 

.90 

59.5 

01 

13 

1 

27.50 

56.60 

138.8 

.43 

94.7 

01 

21 

8 

32.20 

51.00 

164.0 

.91 

24.0 

01 

26 

5 

' 33.88 

46.34 

185.5 

.99 

25.4 

01 

31 

5 

34.00 

42.56 

207.9 

.76 

29.8 

01 

36 

5 

32.80 

38.00 

218.5 

.94 

23.8 

01 

48 

12 

25.61 

30.82 

263.3 

.85 

38.3 

01 

49 

1 

25.46 

30.24 

257.7 

.60 

63.2 

02 

00 

11 

15.20 

30.34 

302.2 

.93 

32.7 

02 

00 

0 

15.02 

30.13 

304.7 

— 

84.1 

02 

08 

8 

9.44 

34.00 

331.2 

.85 

26.5 

02 

12 

4 

7.38 

36.42 

346.4 

.80 

26.8 

02 

15 

3 

6.22 

38.81' 

359.1 

.89 

24.9 

02 

15 

0 

5.96 

39.60 

355.2 

-- 

13.4 

02 

22 

7 

5.62 

45.16 

22.1 

.80 

25.6 

02 

27 

5 

6.00 

50.00 

28.5 

.97 

24.0 

*Rudder  execute  from  strip  chart  at  00:20  = 00:34  video  (14  sec  apparent  lag) 


TABLE  10 

Circle  Test  No.  17,  Speed  Setting  = 8 
Tow  Position  from  Video  Tapes 


Position  of 
Center  of  Gravity 

Time  At  (ft) UCG 


(Min: sec) 

(Sec) 

YP 

Xo 

*# 

ft/sec 

8° 

39:42 

39:48 

6 

5.97 

4.00 

30.32 

36.62 

-3.0 

0.0 

1.07 

-17.9 

39:50 

2 

4.00 

35.02 

0.0 

.80 

0.0 

39:52 

2 

4.00 

40.00 

0.0 

2.49 

0.0 

39:57* 

5 

4.00 

39.20 

15.9 

.16 

15.9 

40:04 

7 

4.22 

42.90 

22.9 

.53 

19.5 

40:09 

5 

5.18 

45.40 

37.1 

.54 

16.1 

40:16 

7 

6.30 

49.00 

43.6 

.54 

26.3 

40:20 

4 

7.66 

50.84 

57.1 

.34 

-26.2' 

40:29 

9 

14.72 

53.55’ 

67.0 

.84 

- 2.0 

40:32 

3 

12.14 

54.50 

79.1 

.92 

9.3 

40:34 

2 

19.04 

54.60 

114.7 

3.45 

25.6- 

40:50 

16 

20.23 

56.72 

103.7 

.15 

74.4 

41:08 

18 

27.72 

54.67 

137.3 

.43 

32.0- 

41:08 

0 

27.98 

54.50 

139.2 

— 

16.01 

41:15 

■ 

30.20 

52.22 

146.6 

.46 

10.9 

41:19 

31.06 

51.14 

159.4 

.35 

17.9 

41:32 

■ 

34.24 

46.16 

179.6> 

.46 

32.1 

41:48 

34.64 

38.31 

203.7 

.49 

26.6 

41:48 

0 

33.94 

38.73 

204.8 

-- 

83.8.' 

42:08 

20 

29.54 

29.48 

232.7 

.51 

27.3 

42  :09 

1 

29.60 

30.78 

245.4 

1.30 

62.8 

42:16 

7 

26.38 

29.36 

256.7.. 

. :>0 

10.3 

42  :19 

3 

25.07 

27.00 

270.0 

.90 

61.0 

42:32 

13 

18.70 

28.24 

289.5 

.50 

8.5 

42:44 

12 

13.45 

29.76 

314.8 

.46 

28.7 

42  :S1 

7 

10.74’ 

31.68 

326.8 

.47 

21.5 

wznm 

3 

9.45 

32.61 

332.3 

.53 

26.5 

2 

9.08 

45.18 

332.1 

6.29 

-26.2 

EI'l 

5 

7.65 

35.02 

338.6 

.60 

15.3 

43  :02 

1 

7.31 

37.55 

339.6 

2.55 

12.7 

EE1 

2 

6.75’ 

36.67 

349.4 

.52 

22.2 

2 

6.29 

37.34’ 

354.5 

.41 

29.0 

wmm 

5 

5.70- 

39.68 

1.5 

.48 

15.6 

jfltl 

0 

5.52 

39.76 

.7 

-- 

1.  £1 

11 

5.80 

45.00 

30.6 

.48 

27.5 

•Rudder  execute  from  strip  chart  at  39.34  = 39.57  (23  sec  apparent  lag) 


TABLE  11 


Comparative  Table  of  Linear  Hydrodynamic  Coefficients 
and  Computed  Rates  of  Turn 
(coefficients  xlO"^) 


Symbol 

Free  Running 

Model  8-ft  Tow 
(1:701 

Rotating  Arm  Tests  of 
Coastal  Barges  (3B+T) 
Models  (1:50) 

PMM  Towing  Tank  Tests 
of  Typical  River  Tow 
(1:12) 

N' 

V 

+ 0.096 

-0.5815 

-0.185 

N' 

r 

- 0.59 

-1.099 

-0.328 

N<5 

- 0.98 

-0.6215 

-0.100 

Y' 

V 

- 15.53 

-3.329 

-1.580 

Y’ 

r 

+ 0.096 

+0.0317 

+0.619 

Y6 

+ 1.63 

+0.7444 

+0.198 

m' 

3.564 

1.1567 

1.457 

L 

7.967  ft 

687.19  ft 

745.06  ft 

6 

O 

O 

CM 

- 20° 

O 

O 

CM 

1 

R/L 

1.81  (2.35)* 

3.44 

5.35 

r ' 

0.554  (0.425)* 

0.291 

0.187 

•Computed  using  y data  from  physical  model  results 


APPENDIX  A:  CALCULATION  OF  HYDRODYNAMIC  COEFFICIENTS 


1.  Mathematical  models  for  ship  motion  calculations  require  the  speci- 
fication of  several  hydrodynamic  coefficients.  Generally,  these  coef- 
ficients have  been  obtained  by  captive  model  testing  in  towing  tanks  and 
theoretical  computations.  The  planar  motion  mechanism  (PMM)  , rotating 
arm  and  other  methods  for  measuring  these  coefficients  are  explained  in 
books  on  naval  architecture . ^ The  review  of  available  information 
revealed  that  maneuverability  tests  had  not  been  conducted  for  tow  flo- 
tillas. Discussions  and  correspondence  with  some  of  the  leading  towing 
tank  facilities  and  towboat  manufacturers  did  not  produce  hydrodynamic 
coefficient  or  tow  maneuverability  data. 

2.  During  the  study,  the  mathematical  model  of  tow  maneuverability  by 

2 

Petrie  became  available.  A method  of  estimating  the  linear  hydrodynamic 
coefficients  was  presented  by  Petrie.  Higher  order  (nonlinear)  terms 
were  ignored  in  his  model,  probably  due  to  the  unavailability  of  mea- 
sured tow  flotilla  forces  and  moments.  The  computation  method  presented 
by  Petrie  was  used  for  calculating  the  hydrodynamic  coefficients  for 
this  study. 

3.  The  following  equations  were  adapted  from  Petrie,  based  on  a linear 
mathematical  model  (see  eqs.  3 and  4 in  the  main  text  of  the  report) : 

(1)  Yaw  moment  equation 

i ii  ii  iii  i 

N = N v + N r + (N.-I  ) r + Nx6  (Al) 

v r r z o 

(2)  Sway  force  equation 

t it  it  i tt  i 

Y = Y v + Y r + (Y.  - m ) t + Y.6  (A2) 

v r v 6 

(3)  Axial  or  surge  force  equation 

x'  = x'  + x'fi  + x'  (A3) 

oo  p 

Terms  not  listed  in  these  equations  compared  to  equations  3 and  4 in 

2 

the  main  text  were  assumed  to  be  zero  by  Petrie.  The  formulas  derived 


Al 


by  Petrie  to  calculate  these  hydrodynamic  coefficients  arc  listed 
below : 

Nv  CdbHbLb  (CCaft  * 7 V - CdtHtLt  tCGfwd  - 7 V 

Nv ' i sr.  " — 


\ PL  U 


(i^  ♦ Lty 


(A4) 


1 ,4,. 

2 Pi  U 


cdb"b  «cc,ft>Lb>3  - CG,ftS1/5 

'h,  * Lt)4 


(A5J 


v 


| pl2u 


“WW  ♦ 

<4,  * Lt>2 


(A6) 


(Nf  - V Cab\  l <“.ft  * V3  ' CGaft3!'3 

N . - I = = - = 

^ 1 r r /■»  . r •»  J 


1 T 5 

2 PL 


Lb  < W 


CatVt  (CCaft5  * CGf.d3^5 


(A?) 


„•  ' (Y*  - ">  KcabV  * <CatVl 

Y . - m = = - 

V i t n . t \ ^ 


1 T3 

2 PL 


+ Lt>' 


(A8) 


i i 

Y = N 
r v 


(A9) 


The  following  values  were  used  in  these  calculations: 

Cdt  - l-S 


Cdb  = °-5 
Cat  " °*45 


Cab  - °-45 

The  resistance  and  propeller  thrust  forces  were  obtained  from  still 


A2 


water,  steady  ahead  tow  tests,  assuming  steady  tow  resistance  is  equal 
to  propeller  thrust.  The  resistance  and  thrust  coefficients  are  given 
by: 

. F . 

(A10) 


X = T 
o 1 


2 2 
2 pLU 


X = 


P 1 T2,!2 

2 pL  U 


(All) 


Rudder  coefficients  were  obtained  from  the  theoretical  lifting  wing 


relations  presented  on  page  518  of  Mandel's  chapter  in  Comstock's  book. 
The  formulas  used  were: 


1 


YS  " 36 


3Y6  AR(1r/2) 


(A3  2) 


V 


N'  . ^aft ; 

6 (Lt  + V 


(A13) 


i • 

X = Y 
6 5 


(A14) 


The  dimensionless  forms  of  the  tow  mass  and  moment  of  inertia  are  given 
below: 


I = 


z 1.5 
2 PL 


* {T  IKb  * <CGaft  * h,  - CGb>2!  * 


t [Kt2  + (CGt  - CGaft)2)}/  } p(Lt  + Lb)5  (A15) 


g 


m = 


(Ab  + ^t^/g 


1.3  1 ,.  . .3 

2 pL  2 p(Lt  V 


(A16) 


4.  The  results  of  these  calculations  are  presented  as  Table  11.  Several 

. 2 


formulas  presented  in  the  text  by  Petrie  were  corrected  based  on  a 
comparison  with  his  computer  code.  It  is  believed  that  the  equations 


A3 


for  the  sum  of  moment  of  inertia  plus  added  moment  of  inertia  and  the 
mass  and  added  mass  are  in  error.  It  is  known  that  in  general  these 
values  should  be  about  twice  the  tow  moment  of  inertia  and  tow  mass. 
Table  1 shows  two  alternate  ways  for  correcting  for  this  apparent  error. 
In  addition,  the  rudder  effects  were  adjusted  to  give  better  agreement 
to  the  measured  model  turning-circle  data.  This  is  further  discussed 
in  the  body  of  the  report . 

2 

5.  Petrie  presented  a method  of  accounting  for  the  increased  velocity 
past  the  rudders  located  in  the  propeller  race.  Calculations  using 
these  formulas  were  tried  but  resulted  in  very  large  coefficients  for 
the  free  running  models.  A study  of  the  Mandel*  chapter  indicated  that 
free  running  model  data  gave  better  correlation  if  the  ship  velocity  is 
used  rather  than  the  propeller  race  velocity.  The  added  resistance  of 
the  small  scale  free  running  models  requires  the  propeller  to  provide 
added  thrust  to  achieve  proper  speed.  On  the  other  hand,  the  approach 
flow  to  the  propeller  is  influenced  to  a large  extent  by  flow  separa- 
tion from  the  model  due  to  added  resistance.  The  two  effects  are 
believed  to  be  compensating  to  some  degree.  The  factor  of  1.72  used  to 
achieve  better  results  with  the  observed  model  turning  tests  would 
indicate  the  propeller  would  have  an  effect  on  the  rudder  coefficients 
of  about  30  percent  (/l .72  = 1.31)  greater  than  the  tow  speed. 


APPENDIX  B:  DEFINITION  OF  SYMBOLS 


m 

v 

v 

r 

r 


U 

* 

6 


NOTE 

P 

g 

Y 

Cb 

^ n 
'at’  ab 


Cdt,Cdb 


CG 

CG 


fwd 


aft 
L 

VS 


mass  of  the  tow 

tow  acceleration  in  the  sway  direction 

component  of  tow  speed  in  the  sway  (y-axis)  direction 

tow  rotary  acceleration  about  the  vertical  axis  5$ 

tow  rate  of  turn  about  the  vertical  axis  =1 

cross-sectional  area  of  rudder  in  sq  ft 

magnitude  of  flotilla  velocity  vector 

yaw  or  heading  angle,  relative  to  fixed  axis,  + from  axis 
to  heading 

drift  angle  or  orientation  of  heading,  relative  to  velocity 
vector,  'Z  from  velocity  to  heading  (B=-sin  yy) 

rudder  angle,  of  the  tow 

thrust  coefficient  per  tow  screw 

total  tow  resistance/2 

2 2 

1/2  pLU 

volume  of  displacement  of  towboat  in  cu  ft 

volume  of  displacement  of  tow  flotilla  in  cu  ft 

displacement  (or  weight  in  air)  of  tow  flotilla  in  lbs 
(or  prototype  tons) 

displacement  (or  weight  in  air)  of  tow  flotilla  in  lbs  (or 
prototype  tons) 

A=yV 

2 1 

density  of  water  in  lb-sec  /ft 

2 

acceleration  of  gravity,  ft/sec 

3 

unit  weight  of  water,  lb /ft 
block  coefficient 

added  mass  coefficient  of  tow  flotilla  and  towboat, 
respectively 

damping  coefficient  of  tow  flotilla  and  towboat,  respectively 

location  of  the  tow  center  of  gravity  from  the  bow  of  the  tow 

location  of  the  tow  center  of  gravity  from  the  towboat/barge 
interface 

total  length  (L  +L.  ),  in  ft 

length  of  tow  flotilla  and  towboat,  respectively,  in  ft 


APPENDIX  B (Cont.) 


VBb 

vs 


draft  of  tow  flotilla  and  towboat,  respectively,  in  ft 

beam  of  tow  flotilla  and  towboat,  respectively,  in  ft 

the  radius  of  gyration  of  the  tow  flotilla  and  towboat,  respec- 
tively, in  ft 

time  constant  of  yaw  acceleration  term  for  ship  or  tow  in  yaw 
equation 

time  constant  of  drift  acceleration  term  for  ship  or  tow  in 
drift  equation 

rudder  constant  of  ship  or  tow  in  yaw  equation 

rudder  constant  of  ship  or  tow  in  drift  equation 

resistance  force 
propeller  thrust  force 

axial  or  surge  force  on  the  tow  in  the  x-axis  direction 

hydrodynamic  yaw  moment 

moment  of  inertia  referred  to  z-axis 

component  of  tow  speed  in  the  axial  or  x-axis  direction 

derivative  of  hydrodynamic  yaw  moment  with  respect  to  yaw  rate 

derivative  of  hydrodynamic  yaw  moment  with  respect  to  yaw 

acceleration 

derivative  of  hydrodynamic  yaw  moment  with  respect  to  sideslip 
velocity 

derivative  of  hydrodynamic  yaw  moment  with  respect  to  sway 
acceleration 

derivative  of  hydrodynamic  yaw  moment  with  respect  to  rudder 
angle 

radius  of  turning  circle  test,  ft 

hydrodynamic  lateral  or  sway  force  component  on  the  tow 
in  y-axis  direction 

derivative  of  hydrodynamic  force  component  in  y-axis  direction, 
with  respect  to  yaw  rate 

derivative  of  hydrodynamic  force  coefficient  with  respect  to 
yaw  acceleration 

derivative  of  hydrodynamic  force  component  in  y-axis  direction, 
with  respect  to  sideslip  velocity 

derivative  of  hydrodynamic  force  component  in  y-axis  direction, 
with  respect  to  sideslip  acceleration 
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Y^  derivative  of  hydrodynamic  force  component  in  y-axis  direction, 
with  respect  to  rudder  angle 


Dimensionless  Forms 

The  dimensionless  form  of  a quantity  is  indicated  by  the  prime  of  that 


quantity.  Examples  are  shown  in  the  following: 


Typical 

Typical  Dimension- 


Quantity 

Symbol 

less  Form 

Force  - x-axis  direction 

X 

X' 

= 

X / 1 l2u2 

- y-axis  direction 

Y 

X' 

= 

Y / I LV 

Moment 

N 

N' 

= 

H / | L3U2 

Mass 

m 

m' 

= 

m/|  l3 

Angular  velocity 

r 

r' 

= 

rL 

U 

Static  force  rate 

Y 

V 

Y' 

V 

= 

Y /|  L2U 
v/  2 

Static  moment  rate 

N 

N' 
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